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Abstract 

The sensitivity of the Statistical Multifragmentation Model to the underlying statistical assump¬ 
tions is investigated. We concentrate on its micro-canonical, canonical, and isobaric formulations. 
As far as average values are concerned, our results reveal that all the ensembles make very similar 
predictions, as long as the relevant macroscopic variables (such as temperature, excitation energy 
and breakup volume) are the same in all statistical ensembles. It also turns out that the multiplic¬ 
ity dependence of the breakup volume in the micro-canonical version of the model mimics a system 
at (approximately) constant pressure, at least in the plateau region of the caloric curve. However, 
in contrast to average values, our results suggest that the distributions of physical observables are 
quite sensitive to the statistical assumptions. This finding may help deciding which hypothesis 
corresponds to the best picture for the freeze-out stage. 

PACS numbers: 25.70.Pq, 24.60.-k 
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I. INTRODUCTION 


The thermodynamics of the multifragment emission in nuclear reactions has 
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The controversy related to the determi¬ 


nation of this quantity is mainly due to the great difficulties involved in the measurement of 


the excitation energy and the breakup temperature of the fragmenting system [23- There¬ 
fore, the existence of a liquid-gas phase transition in nuclear matter could not be clearly 
established or ruled out yet. Furthermore, several statistical calculations, supported by ex¬ 


perimental analyses suggest that the heat ca pac ity at constant pressure. Up, can assume 
negative values under certain conditions |la lla |20|, |3^ whereas other considerations 
have led to opposite conclusions QQ. 

Theoretical aspects related to the smallness of nuclear systems have also been discussed 
and could influence conclusions drawn from statistical calculations based upon different 
constraints y. More specifically, the equivalence of statistical ensembles is expected to be 
valid only for very large systems, i.e. in the thermodynamical limit. Of course, this condition 
is not fulfilled in the reactions we discuss here. In spite of this, statistical approaches are very 
successul in quantitatively describing many features of nuclear multifragmentation |2fll. l21|. 


In this work, we compare the predictions of the Statistical Multifragmentation Model 
(SMM) 


H y 0. co„ce„t.at.g o. .ts n.c.o-ca.o„ica.. caaon.ca. and .sobanc ve..ons. 
By comparing them, one actually confronts two different pictures for the freeze-out stage. In 
the canonical and in the micro-canonical frameworks, the break-up volume is fixed, so that its 
properties are dictated only the by range of the nuclear forces. The underlying assumption 
in the isobaric ensemble is that the freeze-out stage is associated with a given collision rate 
and, therefore, the break-up volume is allowed to fluctuate from one microstate to another. 
Since the freeze-out configuration can be rather different in both cases, one expects to find 
characteristic fingerprints for some observables. 

In sect. m we give a detailed review of the different versions of the SMM used in this 
work and, in particular, of the isobaric ensemble. The comparison among their predictions 
is discussed in sect. uni We conclude in sect. na with a brief summary and with suggestions 
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for experimental analysis. 


II. THE MODEL 


The SMM is based on a scenario in which a hot and compressed source expands, un- 
dergoinga prompt statistical breakup at a density smaller than that of normal nuclear 


matter 


36|. The different fragmentation modes {/} are weighed accorc 


ing to their statis- 
, Isa bs] , although 


tical factors. The latter depend on the statistical ensemble adopted j21 
the essential ingredients of the model do not vary. We present below the three statistical 
versions used in this work. 


The Microcanonical Ensemble 


In the microcanonical version 




36, 


37| , each fragmentation mode / of an excited source 


with mass and atomic numbers Aq and Zq, respectively, must strictly obey the following 
mass, charge and energy constraints: 


^0 — ^ ) 
{AZ} 

Zo = ^ n{^Z , 

{AZ} 

and 


( 1 ) 

( 2 ) 


-Ba,,z, + E* = (w) ^ + E (3) 

^0 ^ ^ {A,Z} 

where ^ represents, within a fragmentation mode /, the multiplicity of a fragment with 
mass and atomic numbers A and Z, respectively, Bao,Zo stands for the binding energy of the 
source, E* is the excitation energy deposited into the system and T its breakup temperature. 
The first term on the right hand side of the equation above corresponds to the Coulomb 
energy of a homogeneously charged sphere at breakup volume V, whereas Vq is the source’s 
volume at normal nuclear density. The fragment energy Ea,z{T, V): 


EaAT, V) 


-Ba,z + ElAT) + exat) + ElAV) 


( 4 ) 
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has contributions from its binding energy Ba,z, translational motion internal exci¬ 

tation energy E\^{T), and from the remaining Coulomb terms, E^^{V), which contribute 
to the total Coulomb energy of the system in the Wigner-Seitz approximation 


Throughout this paper, we use the Liquid Drop formula adopted in ref. 


36, 
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for all fragments with A > 4. Empirical binding energies are used for lighter nuclei. The 
parameters entering in the expression above are Wq = 16.0 MeV, /3o = 18.0 MeV, Cc = 
0.737 MeV, TLasym = 30.0 MeV, and Qasym = 35.0 MeV. The remaining energy terms in Eq. 
(JD) read: 
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izm = , 
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and 
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where eo = 16.0 MeV and = 18.0 MeV is the critical temperature above which the surface 
tension vanishes. Light nuclei, A < 5, are treated as particles without internal degrees of 
freedom. Therefore, they only contribute to the total energy through the binding and kinetic 
energies, except for alpha particles in which case one retains the bulk term in Eq. ((7j) to 
account for its particle-stable excited states. 

The energy conservation stated by Eq. dSl) allows one to determine the microcanonical 
temperature Tf for each fragmentation mode /. Therefore, Tf fluctuates from one partition 
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to the other. The statistical weight tc/ of a fragmentation mode is given by its corresponding 
number of states: 


Wf = exp {Sf) , (10) 

where Sf stands for the entropy, which is obtained from the sum of the contributions due to 
each fragment. It is related to the total energy E and the Helmholtz free energy F through 
the standard thermodynamical expression: 

F = E-TS (11) 

and 

S(T,V)^~(T,V). (12) 

The Helmholtz free energy for a fragmentation mode / can be written as: 


where 


Ff{T, V) = F}{T) + V) + Ff{V) 
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In the above eqnations, gA,z stands for the the spin degeneracy factor of the fragment. It is 
assnmed to be eqnal to nnity for fragments which have internal degrees of freedom whereas 
empirical gronnd state valnes are nsed for the others. The free volnme is related to V by: 


Vfree = V" - Vq = VqX, 0 < X < OO (18) 

and At = \/27rh?/mnAT, where rUn corresponds to the mass of the nncleon. 

The mean valne of any physical observable (O) is obtained by weighing its valne in each 
partition by the corresponding nnmber of states Wf\ 


(O) 


Micro — 


Y.fWfOf 

EfWf 


(19) 


In particnlar, the average pressnre (P) can be caicniated throngh: 


Pf = 


dPf 


(M, - 1)T 0 / v/'’ 


dV V,,.. ' 3 V'l/f ■ 

The kinetic contribntion is proportional to Mf — 1, rather than to the total mnltiplicity, 
Mf = 'Y1,{az}^a,Zi as the center of mass is kept at rest. Even if it is not explicitly stated 
in the formnlae, this condition is consistently imposed in all the versions of the model. 
Therefore, Eq. ©, for instance, is actnally modified to take it into acconnt. 


Since, for large nnclei, the nnmber of different fragmentation modes is hnge [37|, we make 
a Monte Carlo sample of the most important partitions. As the sampling method does not 
select partitions with eqnal probability, the statistical weight w/ is corrected for this bias 
j^|. This procednre is also adopted in the canonical and in the isobaric ensembles presented 
below. 

Finally, a mnltiplicity dependent free volnme has been introdnced in refs. |^. I^|: 


Xf = 


D 




'0^0 


- 1 


( 21 ) 


where D = 2.3ro. This canses the breaknp volnme to flnctnate from one partition to the 
other. 
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B. The Canonical Ensemble 


If the breakup temperature of the system T is hxed, instead of the excitation energy of 
the source, besides its breakup volume, mass and atomic numbers, the canonical ensemble 
is the best suited statistical treatment. The partition function associated with a given 
fragmentation mode reads: 


Zf(r.l/) = exp 


-Ff(T,V) 

T 


Then, the average value of a physical observable is calculated through: 


( 22 ) 


Y.,z?o, 

(0)c = ^ ' ' . (23) 

Since the same essential ingredients are used in all the statistical ensembles, the Helmholtz 
free energy F is also given by Eq. (HD- Therefore, Eq. (EHD still holds for the pressure in a 
given fragmentation mode. The total energy of the system can be calculated though Eqs. 
m and so that it is given by: 


Ej = ^T{M, - 1 ) + ^ 1-Baz + E\z{T)] 

{AZ} 

(24) 

which now fluctuates from one fragmentation mode to the other since the temperature is 
hxed. The excitation energy is then given by: 

E*f = Ef + BaoZo ■ (25) 


C. The Isobaric Ensemble 

If one now Exes the pressure P, instead of the breakup volume, the canonical ensemble 
can be modihed to keep the pressure constant for all microstates. In this case, the partition 
function becomes: 



dV Zf{T, V) exp 



(26) 
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which, by using Eqs. ()lHlllb|l . can be rewritten as: 


Zf%T,P) = Zf{T,P)If{T,P). 


In the equation above, we have dehned: 
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Similarly to the former case, average values are obtained through: 


( 27 ) 


(28) 

(29) 

(30) 


(31) 


E/ zfo, 
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(32) 


The appropriate thermodynamical function for the isobaric ensemble is the Gibbs free 
energy G: 


G = F + PV , 

which is related to the partition function by: 


(33) 


G,{T,P) = -Tlog[Zf°{T,P)] . (34) 

From it, one may obtain the average volume of the system, for a given fragmentation mode: 
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If(T,P) • 
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The energy of the system corresponding to the breakup channel / can then be written as: 


so that 
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(39) 


Owing to the Coulomb factor (1 + x)~^^^i all the integrals above have to be evaluated 
numerically for each fragmentation mode. Since the integrands are bell shaped, the gaussian 
quadrature method is very efficient. 

Finally, from the above expressions, the entropy can be evalutated from the relation: 


O Ef + PV’ - Gf 

^f - -T- 


(40) 


III. COMPARISON AMONG THE ENSEMBLES 

All the example calculations in this work are carried out for a A = 168 and Z = 75 system. 
The same breakup temperature is used in the isobaric and in the canonical calculations. 
We hxed the pressure, in the isobaric ensemble, at p = 0.114 MeV/fm^. Although our 
conclusions are not qualitatively affected by this particular choice, the results below show 
that this is a reasonable value and lies within the expected range for this system Q . Since 
the micro-canonical ensemble requires the energy as input, we use the average excitation 
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energy obtained in the isobaric ensemble, calculated through Eqs. (jdTHdSj) and for a 
given temperature. The breakup volume, which is a free parameter for the canonical and 
micro-canonical ensembles, is chosen in different ways, as explained below. 



E* / A (MeV) 


FIG. 1: Caloric curve predicted by the different ensembles. Upper panel: the breakup volume is 
kept constant for all energies for the canonical and the standard micro-canonical ensembles. Lower 
panel: the average breakup volume predicted by the isobaric calculation is used as input for the 
other ensembles. For details, see text. 


Before we present the model results for the different ensembles, it is important to stress 
that the error bars in all the hgures correspond to the standard deviations of the distribu¬ 
tions, rather than to statistical fluctuations. Indeed, all the curves show a smooth behavior, 
even when the error bars are appreciably large. We included this information in all the 
pictures so as to provide more information about the distributions instead of only mean 
values. 

We start by analysing the caloric curve predicted by the isobaric ensemble, which is 
depicted by the circles in Fig. [T] It shows that, at low temperatures (T < 6 MeV), the 


caloric curve follows the standard Fermi gas expression 


represented by the dashed line in 


this picture), in agreement with previous calculations 3^. However, a sudden change takes 
place at T = 6 MeV and one observes a fairly wide plateau, which signals a liquid-gas 
phase transition. For higher temperatures, the system tends to behave as a free gas, as 


also observed in ref. 


We have checked that the critical temperature depends on the 
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pressure, but the essential features of the process remain the same. 

It is well known that SMM predicts the existence of a plateau in the caloric curve only 
if the breakup volume, for a hxed excitation energy, is allowed to vary from one fragmen¬ 


tation mode to the other j3J, |d 


41 


42j |. Indeed, we also show, in the upper panel of 


Fig. □ the results obtained using two versions of SMM. The one for which the breakup 
volume is determined, for each fragmentation mode, by Eqs. (HHD and (EH), is labelled 
“M.D. micro-canonical” (Multiplicity Dependent), to distinguish it from the “standard” 
micro-canonical version (with no additional labels) in which the breakup volume is hxed for 
all partitions at a given excitation energy. The results clearly show that the temperature 
increases monotonously if the breakup volume is kept hxed for all energies (diamonds) al¬ 
though a behavior similar to that predicted by the isobaric ensemble is observed in the other 
case (squares). In all calculations in which the breakup volume is the same for all energies, 
we adopt D/Vq = 3, for both micro-canonical and canonical ensembles. 

The predictions made by the canonical version of SMM are represented by the triangles 
in the upper panel of Fig. ^ An excelent agreement between the canonical and the mi- 
crocanonical ensembles is obtained, showing the consistency of the model, as long as the 
relevent macroscopic quantites, such as breakup volume, average temperature/excitation 
energy, and pressure, are similar. 

This expectation is reinforced through the comparison among the caloric curves predicted 
by the isobaric, standard micro-canonical, and canonical ensembles shown in the lower panel 
of Fig. m In this case, the breakup volume for the latter two ensembles, for each tempera¬ 
ture/excitation energy, is taken to be the average value obtained in the isobaric calculation. 
It is clear that all the curves collapse into a single one. This implies that the most appropri¬ 
ate statistical scenario for nuclear multifragmentation may not be singled out by examining 
average values only. 

To further clarify the aspects discussed above, we now investigate the average pressure 
predicted by the different versions of the model. The results shown in Fig. |21are calculated 
from Eqs. EOD-dH and Eqs. EnD-E3D, for the micro-canonical and canonical ensembles, 
respectively. A steady increase of the pressure as a function of the excitation energy is 
observed, in the upper panel of this hgure, for the standard micro-canonical and the canonical 
calculations. Since the breakup volume is the same for all excitation energies, this behavior 
is quite reasonable on physical grounds. On the other hand, fluctuations of the breakup 
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FIG. 2: Average pressure predicted by the different ensembles as a function of the excitation energy. 
In the upper panel the breakup volume is kept constant for all energies, except for the M.D. micro- 
canonical version. The results displayed in the lower panel are obtained by using the average 
breakup volume predicted by the isobaric ensemble as input to the canonical and micro-canonical 
calculations. 

volume, from one fragmentation mode to the other, introduced by Eqs. m and (EH), 
lead to a weak energy dependence in the plateau region, as revealed by the results obtained 
employing the corresponding micro-canonical model, represented by the squares in the upper 
panel of this picture. 

The role played by the volume change is emphasized by the results shown in the lower 
panel of the same picture. In this case, the breakup volume for the canonical and micro- 
canonical calculations is chosen so that it corresponds to the average value obtained in 
the isobaric ensemble for each temperature. The pressure now stays fairly constant for 
E*/A > 5 MeV. For lower excitation energies, it decreases owing to the influence of the 
large remnant, which, in this energy range, is almost always present. 

The prediction of a plateau in the caloric curve and the similarity of the average pressure, 
in this region, found in the isobaric and the M.D. micro-canonical calculations, seem to 
indicate that the ad hoc multiplicity dependence, in the latter, Eqs. m and EH), leads to 
similar breakup volumes in both statistical ensembles. This is illustrated in Fig. IHl which 
displays the average breakup volume obtained with these two versions of the model. As 
expected, the predictions are very similar, although the micro-canonical ensemble gives 
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smaller average breakup volumes at high energies. This explains why the average pressure 
increases for E* jA > 10 MeV. A striking feature of the isobaric ensemble is the linear 
energy dependence of the breakup volume after the onset of multifragment emission. This 
is expected to happen for a gas with no internal degrees of freedom but, particularly at the 
plateau, the contribution from complex fragments is not negligible. 
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Isobaiic 

M.D. Micro-canonical 


2 - 
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10 


E / A (MeV) 




15 


FIG. 3: Average breakup volume as a function of the excitation energy predicted by the isobaric 
and the micro-canonical ensembles. The dashed line just emphasizes the fairly linear dependence 
of the breakup volume on the excitation energy in the isobaric ensemble. For details see text. 


Although the compatibility among the different versions of the model can be understood, 
some qualitative aspects still depend strongly upon the version one uses. This is demon¬ 
strated by the plot of the entropy as a function of the temperature in Fig. 0] In the upper 
panel, the canonical and the standard microcanonical calculations were carried out at hxed 
breakup volume for all temperatures. One sees that the entropy rises smoothly as the ex¬ 
citation energy increases. However, a backbending is observed in the results obtained with 
the M.D. microcanonical version. This led the authors of ref. Q to predict negative heat 
capacities in the plateau region. However, our results (Fig. |2I) show that the pressure is 
only approximately constant in the microcanonical calculations. Thus, one should be care¬ 
ful when drawing conclusions about Cp from the entropy plot within this microcanonical 
calculation, since the microstates are not subject to a constant pressure. 

On the other hand, from the results obtained with the isobaric ensemble, depicted by the 
circles in this picture, one can certainly assert that Cp is very large, but always positive, 
in the plateau region. Nevertheless, the consistency of the model is not affected by this 
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FIG. 4: Entropy per nucleon as a function of the temperature. Results corresponding to the 
canonical and the standard microcanonical ensembles are obtained at a constant volume for all 
temperatures in the upper panel, whereas the average values given by the isobaric ensemble were 
used as input in the lower panel. For details see text. 


disagreement as is shown in the lower panel of Fig. 01 in which the breakup volume for the 
canonical and the microcanonical calculations are taken from the average values predicted by 
the isobaric ensemble at each temperature. Also in this case, a very good agreement among 
the ensembles is found. The microcanonical calculations still seem to predict negative Cp 
in a narrow temperature region, although one should keep in mind the remarks above. The 
relevant point here is that all ensembles make very similar predictions under equivalent 
conditions. 

Even though it is difficult to determine experimentally the heat capacity, other observables 
are very sensitive to the statistical assumptions and could, therefore, help to decide which 
picture is more convenient to describe the multiframent emission mechanism. 

For instance, in Fig. El we show the isotopic distributions for a few selected cases obtained 
using the different versions of the model. One notices that the isobaric, canonical, and the 
standard microcanonical ensembles give very similar isotopic distributions. As above, the 
breakup volume for the canonical and the microcanonical models is taken from the isobaric 
calculations. In contrast, the isotopic distributions predicted by the M.D. microcanonical are 
systematically narrower. Since precise data are avaliable for many reactions important 
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FIG. 5: Isotopic distributions predicted by the different ensembles. The yields are normalized 
relative to the maximum for each Z value. 

conclusions may be drawn from the study of this observable. 



FIG. 6: Energy distribution obtained with the isobaric and canonical ensembles calculated at 
T = 6 MeV. For details see text. 

Another key observable is the excitation energy distribution, which is shown in Fig. 
El for the isobaric and canonical ensembles at T = 6 MeV. As can be noticed, it turns 
out to be quite sensitive to the statistical assumptions. The canonical calculations have 
been performed in two different ways. In one case, the breakup volume corresponds to 
the average value obtained in the isobaric calculation and is labelled canonical-V. In the 
second canonical calculation, the breakup volume is chosen so that it gives the average 
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pressure p = 0.114 MeV/fm^, and labelled canonical-P. These two canonical calculations 
agree qualitatively, as both give bell shaped distributions, which essentially differ by the 
mean values. On the other hand, when all the microstates are at the same pressure, the 
energy distribution has a narrow peak at low excitation energy and a smeared bump at high 
excitation energies, with a large gap between them. The low energy part of the distribution 
correspond to the presence of a large remnant (the liquid phase) in the breakup channel, 
whereas essentially light fragments are found at high energies (the gas phase). 

This is illustrated in Fig. [71 which shows the results obtained with the isobaric and 
canonical ensembles. The canonical calculations are done in the manner just described. 
One should notice that the tendency to produce a large remnant in the breakup channel 
seems to be a pressure efffect as the canonical calculation which mimics the isobaric ensemble 
appears to produce heavy fragments and even shows a small bump near Aq. This behavior 
is not seen in the other canonical calculation. 



FIG. 7: Distribution of the largest fragment within each fragmentation mode obtained with the 
isobaric and canonical ensembles at T = 6MeV. 


IV. CONCLUDING REMARKS 

We have examined the sensitivity of the predictions of the SMM to its statistical as¬ 
sumptions. We found that the microcanonical, canonical, and isobaric implementations of 
SMM predict very similar caloric curves, and other (average) physical observables, provided 
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macroscopic variables, such as temperature, excitation energy, and breakup volume, are the 
same in all versions of the model. 

Although this is not the main point of this work, it is worth reporting that we found large 
values of Cp at the plateau region of the caloric curve, but no indication of negative heat 


f. Q. 


capacities in our isobaric calculations, in agreement with ref. UJ . We have concentrated 
on the Aq = 168, Zq = 75 system, but we have also checked that this conclusion still 
lolds for lighter systems, such as the ^‘^^Ag nucleus, in contrast with other calculations 


30H 


M 


35|. More specihcally, we observed negative Cp only in the microcanonical 


calculations, as already reported 


in refs. Q,Q,Q- 


Nevertheless, the predictions made 


by the isobaric ensemble should be more reliable since it is the only statistical ensemble 
in which all the microcospic states are subject to the same pressure, whereas, in the other 
cases, this condition if fulhlled only on the average. We should also mention that we did not 
investigate the sensitivity of this conclusion to the model parameters related to the binding 
and excitation energy of the source. As reported in ref. [l^ these parameters may influence 
the mass region in which negative Cp values are obtained in the statistical calculations. 

We would like to stress the point that the distributions of physical observables, rather than 
their mean values, seem to be fairly sensitive to the statistical assumptions. The relevance of 
this conclusion goes beyond the fact that the equivalence of the statistical ensembles should 


hold only in the thermodynamical limit js]. Indeed, some predictions made upon different 
assumptions agree very well, as long as the most relevant input global parameters do not 
differ appreciably. Strictly specking, one may distinguish among the different calculations. 
However, we belive that experiments cannot be precise enough to single out one of them 
only by examining average values. Therefore, although nuclear systems are small and nuclear 
reactions cannot be hne tuned, the distributions of some physical observables should bear 
clear signatures related to the underlying statistical assumptions. Then, the analysis of such 
distributions may be valuable in establishing an appropriate scenario for the freeze-out stage 
in nuclear multifragmentation. 
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